Noname manuscript No. 

(will be inserted by the editor) 



Immo Sollner • Benjamin Gschosser • 
Patrick Mai • Benedikt Pressl • Zoltan 
Voros • Gregor Weihs 

Testing Born's Rule in Quantum 
Mechanics for Three Mutually 
Exclusive Events 

Received: March 6, 2011 / Accepted: July 26, 2011 

Abstract We present a new experimental approach using a three-path in- 
terferometer and find a tighter empirical upper bound on possible violations 
of Born's Rule. A deviation from Born's rule would result in multi-order 
interference. Among the potential systematic errors that could lead to an ap- 
parent violation we specifically study the nonlinear response of our detectors 
and present ways to calibrate this error in order to obtain an even better 
bound. 
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1 Theory 

1.1 Introduction 

Two path interference is often regarded the characteristic phenomenon of 
any wave theory. It was first investigated by Thomas Young in 1801 when he 
tried to answer the question whether light was made of waves or particles. 
A similar experiment using electrons was repeated by Claus Jonsson in 1961 
to show that the interference of two possible paths exists for electrons, as 
predicted by quantum mechanics. This is still often referred to as one of the 
most beautiful physics experiments of all times [1 j. Born's rule allows us to 
calculate probabilities from quantum mechanical wavefunctions and results 
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in interference terms that always originate from pairings of paths, but never 
triples, etc. Of course, multi-path interference will result in more complex 
interference patterns, but only by virtue of the multiple path-pair interfer- 
ence terms. Only very recently have there been investigations, experimental 
and theoretical, of possible new, higher-order phenomena going beyond the 
standard description of interference, when three or more paths are made to 
interfere [2,3,4 . These investigations have been sparked by the seminal work 
of R. D. Sorkin [5]. 

Out of an idea to implement quantum theory as measure theory on space- 
time Sorkin analyzes interference as a deviation from the classical additivity 
(sum rule) of the probabilities of mutually exclusive events. If one considers 
more than two mutually exclusive events or paths A, B, C, . . . one can define 
a hierarchy of higher-order interference terms. 

If -Pa 

T AB = P AB -PA-PB (1) 
tflc = P ABC ~ PAB ~ PAC ~ P B C + P A + P B + P C 

These three terms are the zeroth, first, and second order interference 
terms. If in a general theory the interference term 1^ is nonzero, this is 
indicative of n th order (or (n + 1) path) interference and the theory is said to 
violate the n th sum rule 1^ = 0. Here Pa is the detection probability of an 
interference experiment with paths A, £?, and C when only path A is open, 
Pab when paths A and B are open, and so forth. 

If the zeroth sum rule is fulfilled, the measurement becomes trivial, and 
we will therefore only concern ourself with the case where lf ] ^ 0. The 
next higher sum rule, = 0, holds for classical particles and is violated 

by quantum mechanics and classical waves. The violation ^ means 
that our system exhibits two-path interference. In classical wave theory and 

(2) 

in quantum mechanics the second sum rule holds, I AB c = 0- ^ n ^ ac ^ this * s 
true for any theory where there is a square law relation between the energy 
(or probability) density and the field amplitude [3 . This square law relation 
is of course Born's rule and is a basic postulate of quantum mechanics. 

The purpose of our experiment is a direct test of Born's rule by looking 
for any violation of the second sum rule. The mutually exclusive events are 
passages through three possible paths of an interferometer as shown in Fig. [I] 
In previous versions [2j[3] these events were approximated by passage of single 
photons through either of three adjacent slits in a diaphragm followed by 
detection in the far field. It has been pointed out, for example in Ref. [7], 
that slits, in particular, if they are closely spaced are not entirely independent. 
This crosstalk will, however generally, vanish very rapidly with the separation 
of the respective paths. Therefore in an interferometer with a macroscopic 
path separation of several ten thousand optical wavelengths, the crosstalk 
caused by opening or closing another optical path will be much much smaller 
than what can be detected. 
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Fig. 1: Generic three-path interferometer composed of two three-way beam- 
splitters (also called 6-ports or "tritters" [6]) and a phase shift associated 
with each of the paths A, £?, and C. Since only the relative phase matters, 
we will set 4>b =0 without loss of generality. 

The experimental quantities that are being measured are the count rates 
R x , where x designates the open interferometer paths. The eight possible 

(2) 

combinations of open and closed paths are then combined as in I ABC to 
form the quantity 

e = Rabc ~ Rab - Rac ~ Rbc ^Ra + Rb + Rc- Ro, (2) 

where Rq is the background rate, taken with all interferometer paths closed. 
This is the experimental equivalent of Eq. [I] and should therefore always be 
zero in quantum mechanics. Since we can not measure the probabilities in 
Eq. [I] directly, the equations are written in terms of count rates, R oc P. 
The last term in e has to be subtracted to eliminate the total contribution 
from any background rate that is constant over all eight combinations, for 
example the detector dark count rate. We verified that the detector dark 
counts are indeed constant in time and independent of the history, i.e. of 
signals applied earlier, within our measurement precision. Beyond that it 
is difficult to check whether the background - e.g. stray light - is indeed 
constant. Any background caused by a closed shutter cancels identically by 
virtue of the definition of e. Likewise, imperfect coherence of the paths has no 
bearing on the result. Background contributions that depend in a non-trivial 
way on the positions of all the shutters would cause a non-zero e and are thus 
sources of systematic errors that can only be eliminated by improvements of 
the experiment. 

The quantity e is further normalized by 

8= \Rab—Ra—Rb+Ro\ + \Rac—Ra—Rc+Ro\ + \Rbc—Rb—Rc+Ro\, (3) 

which is the sum of the absolute values of the expected two-path interference 
terms. This newly defined quantity [3], k = f, makes it more meaningful 
to compare the results from different experimental approaches, because it 
normalizes the hypothetical second order interference by the expected first 
order interference. 

In the following we will discuss a new experimental approach to this mea- 
surement, a three-path interferometer, and will focus on one of the major 
reasons for experimental deviation from the expected null result, the nonlin- 
earity of our detector. Using a simple model for the nonlinear behavior of the 
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Avalanche Photodiode (APD) illuminated by a Poissonian light source, the 
deviation from the null measurement can be predicted, compared to experi- 
mental results, and finally used to correct this systematic error. 



1.2 Detector Nonlinearity 

The main reason for the nonlinearity of our detector is its dead-time. Each 
detection event of an APD is followed by a period of time during which the 
detector is unresponsive. During this time the detector can essentially be 
viewed as being turned off. If the efficiency of the detector during the active 
time periods, can be viewed as being independent of factors such as the the 
average count rate or the total incident rate then it can be ignored. Therefore 
the efficiency was set to one for the following discussion. To quantify this 
effect we want to describe the number of photons detected, i? det , as a function 
of the total rate R of incident photons per unit of time. The simple model 
presented below shows that the number of incident photons lost because their 
arrival time is within the dead-time caused by a previous detection does not 
grow linearly with the total number of incident photons, and therefore leads 
to a nonlinear measurement response of the detector. [8] 

An attenuated laser is used as the photon source in this experiment and 
the photon emission statistics can therefore be described by a Poissonian 
distribution so that the average number of incident photons is simply i?TQ 
In Eq.[4]it is shown that the rate of detected photons per second is equivalent 
to the rate of incident photons minus the number of photons that are expected 
to arrive within the dead time r after each detection. 

R det = R- R det Rr (4) 

Rearranging the equation to solve for R det and R respectively yields. 

R R de ^ 

Rdet = YTr^ ' r= i-r^t (5) 

So far we assumed that a detection can only be triggered when a photon is 
incident on the detector. For this model to be more accurate in the low inten- 
sity regime the effect of detector dark counts and other sources of constant 
background detections have to be taken into account. Consequently we can 
define the detector transfer function / and its inverse / as 

«--/WJW:- 1+ y + \ )r (•) 

fi = ;(^):= ^;'^> . (7) 

This model can now be used to predict an expected deviation from the zero 
measurement solely due to the nonlinear detector response. For this, however, 
it will be important to precisely determine the detector dead time and back- 
ground count rate, and to verify that the model does describe the detector 
response. 

1 This holds for any stationary light source if the time interval is larger than the 
source's correlation times. 
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1.3 Calculating ^ det 

In this section we show how the equations describing the detector nonlinear- 
ity, introduced above, are used to predict expected deviations from k = 0. To 
do this we calculate the incident intensities of the multi-path combinations 
from the measured values of the single-path cases R^, R B et , and R^. Then 
the two and three-path incident rates are 

Rab = Ra + Rb + 2^R a R b cos {<t> A ) (8) 

Rabc = R A + Rb+Rc + 2y/R A R B cos (0 A ) 

+ 2y/R A Rc cos {(j>A - <t>c) + 2^/RbRc cos (-^c), (9) 

where <p A and <pc are the phase shifts applied to paths A and C, respectively 
and path B experiences no phase shift, i.e. serves as the phase reference. 
The incident photon rates are calculated using the inverse detector transfer 
function, i. e. R A = /(i?^ et , Rq), and similarly for Rb and Rc- 

From these inferred incident rates we can predict the rates Rj^ B , R^c, 
R B e ci an d R^bci an d thus ft det , taking into consideration the nonlinear re- 
sponse of the detectors and background detections. Here only the detector 
dead time r, the background count rate Rq, and the measured rates R^, 
R B et , and R^ are used as input parameters. We can then compare these to 
the values that we obtain directly from the data of all eight path combina- 
tions to see if there is any amount of second order interference beyond what 
can be explained by the detector nonlinearity. 

Finally we could deduct the systematic error caused by the nonlinearity 
from the measured n value. Any such correction has to be done with caution 
because in general we may not be able to distinguish the nonlinearity of 
a detector from a nonlinearity of nature, i.e. a violation of Born's rule. A 
correction is only possible if the detector nonlinearity is being calibrated by 
a method that does not use interference, for example the beam combination 
method, which is detailed in the following section. 



2 Experiment 

2.1 Testing detector nonlinearity 

For a precise determination of the detector nonlinearity, i.e. the dead time r, 
we set up a direct measurement using the beam combination method. This 
method is based on the addition of non-interfering beams on the detector [9] . 
The setup in Fig. [2] was used for these measurements. For each individual 
measurement two independent laser beams were adjusted in their intensity 
and we measured four values, the dark counts, the count rate for each beam 
separately open and the total count rate when both beams were open. This 
is repeated for a variety of laser intensities in order to be able to map out the 
behavior of the detector for the desired range with increasing accuracy. The 
intensities are chosen by changing the attenuation in the two independent 
arms using variable attenuators. In order to infer the dead time, several 
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Fig. 2: The experimental setup being used for the characterization of the 
detector nonlinearity using the beam combination method. 



assumptions had to be made, yet, fewer than for any other method we are 
aware of. Most importantly we assumed: 

— That the intensities in the two arms will remain constant on the time 
scale it takes to record the four values, for this we have to examine the 
laser fluctuations and fluctuations on transmission of the attenuators. 

— That the dead time has a constant value for the range of incident inten- 
sities investigated. 

— That the detector efficiency, when it is active, is independent of factors 
like the total incident intensity and the average count rate. 

The value of the dead time can be extracted from these measurements using 
a method similar to that presented in Ref. [10 . A more detailed paper on the 
nonlinearity measurement and parameter estimation is under preparation. 

The avalanche photon detector used in this experiment are specified to 
have a dead time of r < 50 ns, and our measurements lead to a dead time of 
r = 47(2) ns and a dark count rate of 284(4) s" 1 . 

In this nonlinearity measurement we use the incoherent addition of beams 
to measure the detectors response to an increased photon (or energy) flux, 
independent of how these rates are connected to the respective amplitudes. 
Therefore we can use this value to correct the systematic error in k to reveal 
any actual violation of Born's rule. In order to make this more convincing, 
we can always vary the overall rate and then check whether our estimated 
systematic error varies accordingly. Any deviation from the expected must 
then either be another systematic error or a violation of Born's rule. 
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2.2 Three-path setup 

We built a three path interferometer using the zeroth and two first orders 
of a transmission grating as the three independent paths. Two of the paths 
pass through independently controlled phase plates, and all three paths are 
back reflected on themselves by a common mirror. The outgoing beam is 
separated from the incoming beam by a double pass through a quarter wave 
plate (QWP) and a polarizing beam-splitter, see Fig. [3j 

The phase plates are 0.9 mm thick antireflection coated glass plates, 
mounted on motorized rotation stages. To compare the measured data to 
the modeled expectations, the rotation angles of the phase plates have to be 
transformed to the actual phases denoted by <pA and 0c, where n\ and n 2 
stand for the refractive indices for air and glass respectively. 



2tt 

T 



2d 



n\ — n 2 



n 2 - ni cos (6 A - 0' A ) 



(10) 



where n2sin#^ = nisin^. The outgoing beam is spatially filtered with a 
pinhole to minimize the contributions from background caused by scattering 
on any of the surfaces in the setup and to allow us to measure a well-defined 
position in the space given by the two phases. Finally the beam is coupled 
into the detector via a multimode fiber. 

A full measurement cycle consist of a 2-dimensional raster scan, using the 
phase plates shown in Fig.[3j with several measurements of n at every position 
in phase space. One n measurement consists of the eight possible combina- 
tions of three shutters independently open and closed. The sequence in which 
these eight combinations are taken for each ft, are randomized to avoid any 
effects caused by slow drifts of the laser. In addition to the complete phase 
space scan we chose a few particular phase settings, where we performed a 
larger number of ft measurements in order to reduce the statistical error. 
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Fig. 3: The experimental setup with the individual shutters for the three 
paths, and the two phase plates to set the position in phase space. 
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2.3 Results 

First, we measured the three different two-path interference visibilities and 
compared them to the achievable maximum values to check the quality of 
alignment. Here it is important to note that path B, the central path, has 
more power than either of paths A and C therefore we couldn't achieve 
high visibilities for the AB and AC path combinations and neither could we 
achieve 100% visibility for the full three-path interferometer [TT] , Then the 
three path interference pattern is measured and compared to the theoretical 
expectations. 




Fig. 4: Contour plots of the three-path interference pattern in two- 
dimensional phase space, where the values on the axes give the phase set 
by the two phaseplates. The crosses mark the phase space maximum used 
for more extensive n measurements. In (a) we show the theoretical model us- 
ing the independent single-path intensities as input parameters. In (b) are the 
the measurement results at the same single-path intensities, i?^ et = 2.08 kcps, 
i?det = 5 75 k C ps 5 i^det _ ^ gg kcps. The only variable used to enhance agree- 
ment between the two plots is a slight adjustment of the absolute starting 
position of the phase plates entered into the theoretical model. This was less 
then 1° and had to be done since the initial starting position could not be 
determined with high enough accuracy. 



A good agreement between theory and experiment is observed, however 
some long time scale stability issues were found when different scans taken 
over several days were compared. To minimize the influence of phase insta- 
bility in the measurement we chose the three path maximum, which in Fig. [4] 
corresponds to m ax = (0c? 0a) = (0.197T, 1.7tt) (marked by a cross), as the 
ideal position in phase space for more accurate k measurements. In Fig. [5a] it 
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becomes clear that this is one of the positions in phase space with the small- 
est gradient and it is also the position with the largest absolute deviation 
for K det due to detector nonlinearity. The hope is that other nonlinear effects 
would also be the largest at this position. 




Fig. 5: (a) Contour plot of the theoretical prediction of tz det in phase space 
using the nonlinear model describing the detector using the independent 
single-path intensities as input parameters. The white cross marks the po- 
sition of the three-path maximum max . In (b) one can see the relationship 
between the intensity at the three-path maximum, i?^^(0 m ax)? an d the 
corresponding calculated ft det (</> max ) (solid line) as well as the values ft exp 
(circles, for error bars see Table 1) measured at the same position. 



Each k measurement consisted of the average of one thousand runs of 
the eight necessary combinations at one fixed position in phase space. As 
mentioned above the three-path maximum was chosen as the preferred po- 
sition in phase space. The phases corresponding to this position were found 
by first doing a fast scan of the entire two-dimensional space like the one 
shown in Fig. [4bJ These angles were then used to take long n measurements 
at this preferred position. Here we were somewhat limited by the accuracy 
with which our motor allowed us to go to a predetermined position and ended 
up measuring for various count rates close to the phase space maximum as 
show in Table [U 



3 Conclusions 

In conclusion our results show good agreement with our model including the 
detector nonlinearity. Future measurements will investigate different areas of 
Fig. [5b] to see if the general predicted trend holds true. When phase stability 
in the interferometer is further improved it will also be of interest to inves- 
tigate different areas of the phase space in Fig. [5a[ We should also point out 
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35925 


-0.0011 


-0.0015 


0.0029 


111288 


-0.0033 


-0.0050 


0.0018 


260934 


-0.0077 


-0.0065 


0.0019 


451121 


-0.0134 


-0.0142 


0.0010 



Table 1: Comparison of the theoretically predicted deviation from k = 0, 
using the detector nonlinear ity model, ft det , and our actual measurement 
results, ft exp for the position in phase space given by ^ max = (0cs0a) = 
(0.197T, 1.77r), which is the intensity maximum of the three-path interference 

Here, An is the standard error of the 



pattern, R^c^ as snown m Fig. 
mean. 
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that the relationship shown in Fig. |5b| depends not only on the dead time of 
the detector but also on the photon statistics of the light source used. The 
ideal case would be the emission from a single quantum emitter with the ex- 
citation events separated by more than a detector dead time. In this scenario 
no effects due to detector nonlinearity would be seen and the detector dead 
time would only limit the maximum achievable count rate. 

Furthermore, with this experiment we can bound the absolute amount of 
multi-order interference as measured by \n\ to be smaller than 0.0015±0.0029. 
This represents an improvement of about one order of magnitude over the 
best previous limit [3]. For this new limit we have not used any correction 
of systematic errors such as the detector nonlinearity model. With a more 
accurate exploration of the relation between count rate and n we expect to 
be able to reduce the bound by several orders of magnitude using the same 
experimental set-up. 
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